Thermal breakage of a discrete one-dimensional string 
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We study the thermal breakage of a discrete one-dimensional string, with open and fixed ends, in 
the heavily damped regime. Basing our analysis on the multidimensional Kramers escape theory, we 
are able to make analytical predictions on the mean breakage rate, and on the breakage propensity 
with respect to the breakage location on the string. We then support our predictions with numerical 
simulations. 
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Recently, there is much discussion on the possibility of 
exploiting biopolymers as functional materials [3, 0, y, 01 • 
To achieve this goal, the stabilities of such materials have 
to be thoroughly investigated. Furthermore, the facts 
that the biopolymers are necessarily finite and consist 
of discreteparts, such as individual peptides in an amy- 
loid fibril [31, have to be taken into consideration. As a 
step towards this direction, we study here a toy model 
for the breakage of a discrete one-dimensional string un- 
der thermal fluctuations, in both fixed-ended and open- 
ended configurations (c.f. Fig.[T]). This problem has been 
studied previously by numerical simulations ^, 0] and 
theoretically with phenomenological assumptions on the 
effect of friction on the collective modes [1, Q. Multi- 
dimensional Kramers escape theory has also been applied 
to the study of breakage in a one-dimensional ring ■ 
The energy profile for the bonds in the string is usually 
modeled by a quadratic potential at the minimum en- 
ergy region, and by an inverted quadratic potential at 
the breakage point. Here, we employ a simplified model 
where all bonds are assumed to be Hookian up to the 
breakage point. This model has the virtue of rendering 
the theoretical analysis asymptotically exact as temper- 
ature goes to zero. By studying in detail the energy de- 
pendency on the collective modes, we are able to employ 
the multidimensional Kramers escape theory to predict 
the breakage rate and the breakage propensity with re- 
spect to the breakage location. These predictions are 
then verified by numerical simulations. 



A. String with fixed ends 

We consider the dynamics of a one-dimensional string 
modeled as a collection of M masses connected by springs 
with identical spring constant k. We also assume that we 
are in the heavily damped regime, i.e., the inertia terms 
are ignored, which is reasonable for many biopolymers 
in typical experimental conditions [ll|. We assume that 
the beads in the string are initially at the minimal en- 
ergy configuration, i.e., each consecutive pair of beads is 



FIG. 1: (a) A ID mass-spring system with fixed ends, (b) 
A ID mass-spring system with open ends. In both cases, 
the separation between the beads are assumed to be of unit 
length. 
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separated by a unit distance. Denoting the positional de- 
viation of the n-th bead from the initial configuration by 
x„, the equations of motion under thermal perturbation 
are of the form: 



dx K , 

dt C 



g 



(1) 



where C is the damping coefficient and g is a Gaussian 
noise such that 
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and 
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Smn6{t - t') 



0-12 



(3) 



i.e., Ann = 2 and An+i,n = An,n+i = -1- 

As A is symmetric, it is diagonalizable by a set of or- 
thonormal vectors. Let D be the diagonal matrix such 
that VDV^ = A where V is the corresponding orthogo- 
nal matrix, we have 
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where p = V^k and the thermal perturbation term g 
remains the same since V is an orthogonal matrix. The 
diagonal elements of D are 



A., 



and the entries in V are: 



M + 1 



Vns = C sin ■ 



M+1 



(5) 



(6) 



In the above equation, C ~ ^2/ [M + 1) is the normal- 
ization factor so that V^^ ~ 1 for all s (c.f. Appendix 

EH). 

The extension/contraction of the n-th spring (desig- 
nating the spring before the n-th bead) is given by: 



ei = 

or in matrix notation, 



[Vv\m - [Vv>\m~i , 



(7) 
(8) 
(9) 



(10) 



where is a matrix of dimension {M -f 1) x Af such that 

(11) 
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l<n< M 



We are interested in the Mean First Breakage Time 
(MFBT) defined as: 

T = inf(t>0 max{|e„(t)|} > 6) . (12) 

t L n J 

In physical terms, we are to find the average waiting 
time before any of the e„ is extended or contracted by 
an amount b where b < 1 (l6j . The above problem is 
equivalent to a multi-dimensional Kramers escape prob- 
le m I12I Il3| , and we will employ the formalism developed 
in [1J| for our analysis. 

In terms of the normal modes {p}, the total energy 
corresponding to the entire string is 



Uip) 



and the exit boundaries are defined by: 

^ WnsPs = ±b , for 1 < n < A/ 



(13) 



(14) 



Our first task is to find the exit routes with the minimal 
energy. To do so, we employ the Lagrange multiplier 
method to minimize the following quantity with respect 
to p: 

t/(p) + z„ WnsPs ±b^ , for 1 < n < M+1 (15) 



FIG. 2: (Color online) Consider the case of having three beads 
connected by one black and one red spring. The energy profile 
in the p-space is depicted in the surface plot in the lower-left 
figure. Breakages of the black (red) spring correspond to the 
black (red) boundaries on the energy plot. A possible trajec- 
tory that leads to breaking the black spring is schematically 
depicted in white. The energy profile at the upper-right exit 
boundary is depicted in the lower-right figure. The MFBT is 
partly determined by the curvature of the potential energy at 
the exit point (c.f. Eq. 1^ ). 
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where z„ is the corresponding Lagrange multiplier. The 
solution to this minimization problem is that for each n, 
there are two minimizing vectors, p("'+ and p^"^~, of the 
forms: 



such that the corresponding energy is: 



(16) 



(17) 



Since (c.f. Appendix lA 2\i 

s 

the minimal energy is the same for ah n: 

C/(P^"^^) = ^^^^^ , l<n<A./ + l. (19) 

As a result, there are 2M + 2 exit points at the boundary 
that correspond to the same energy. We will denote these 
exit points by Q^*^^: 



)W_.^1^fe^^fc/2l,s6(M + l) 



l<k<2M + 2. 

(20) 
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In Fig. [2l we depict the energy profile and the corre- 
sponding exit points for the case of a three-bead system. 

With the formahsm developed in in the asymp- 
totic limit of K — *■ oo (or equivalently, ksT — > 0), the 
MFBT can be expressed exactly as: 



C\/27rfcBT exp 



2MkBT 



where 



= det 



= det 



92 [/ 



dprdps 



Ha. 



(21) 



(22) 



dprdps 



1 < r, s < Af - 1 (23) 



with Pj.{1 < r < M — 1) being a set of basis that 
are perpendicular to the direction of the fc-th exit 
route, i.e., they are perpendicular to the direction 
(W^[fe/2],i, ■ • ■ , W^[fe/2],M)- In other words, (/jq corresponds 
to the Hessian of the potential energy at the origin, and 
corresponds to the Hessian of the potential energy 
within the hyperplane that has its normal pointing along 
the A:-th exit route. Note also that physically, |kZ)Q^*''^| 
corresponds to the magnitude of the potential energy's 
gradient at the fc-th exit point p^ . 

In Appendix lASK Eq. and Appendix |X5] (for Eqs 



and 



:DQ 



(k) 



we argue that 

{ (A/+l)Kb 
CM 
^/2{M+l)K.b 
CM 



M- 



iM/2 



k = 1,2,2M 
otherwise 



k = 1,2,2A/ 
otherwise . 



l,2Af + 2 
(24) 



1,2M 



(25) 



FIG. 3: (Color online) Breakage events for the fixed-ended 
string. Upper plot: The ratios of the MFBT from simulations 
vs. the MFBT from theory. Each marker represents 1000 
samples. The parameters are: = 10, fcaT = 1 and b = 
0.1. Lower plot: The breakage frequency with respect to the 
breakage location for the case of M = 8 with the combined 
data from the cases of k = 1600, 1800, 2000. 
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break half as frequently as the springs in the middle. 
Fig. [3] shows that the analytical predictions are in good 
agreement with simulations. Physically, the fact that the 
extremal springs break half as frequently may be seen 
from the fact that they are connected to the rigid wall 
on one side and so they are subjected to about half the 
amount of thermal fluctuations. 



Substituting the above quantities to Eq. (El 
at the following prediction for the MFBT 



we arrive 




c 



1) K^I%{M + 1) 



exp 



(M + 1)k&2 
2kBTM 



(27) 

Fig. [3] demonstrated that the convergence of the numer- 
ical results to the analytical predictions as k increases. 

Besides the MFBT, this formalism is also capable of 
predicting the propensity for breakage with respect to 
the location of the spring in the string. According to 
Eq. 5.1 in the probability of breakage at the n-th 
segment is given by: 



Pr(n) 



20^^/VjJQ(")/C| 



1/{2M) 
1/M , 



n=l,M+l 
otherwise . 



(28) 
(29) 



This signifies that the breakage propensity is uniform for 
all springs except for the two extremal springs, which 



B. Strings with open ends 

For an open-ended string as depicted in Fig. [T|d, the 
equations of motions are governed by: 



—r- = X2~ Xi - 1+ gi 

at 

dXn 

—7— = Xn+1 - 2Xn + Xn-1 + gn , l<n<M 

at 



(30) 



dt 



XM-1 - Xm + 1 + (JM 



where Xn{t — Q) — n and g is as defined in Eq. ([2]). We 
now let 



M - 1 



l-^J^E^"' (31) 



i.e., these new coordinates are defined in relation to the 
center of mass of the string. With these transformations, 
the equations of motion become 



dx K 



^x-l-g 



(32) 
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FIG. 4: (Color online) Breakage events for the opened-ended 
string. Upper plot: The ratios of the MFBT from simulations 
vs. the MFBT from theory. Each marker represents 1000 
samples. The parameters are: C, = 10, fcaT = 1 and b — 
0.1. Lower plot: The breakage frequency with respect to the 
breakage location for the case of M = 8 with the combined 
data from the cases of «: = 2000, 2200, 2400. 



>. 




o 

0) 


1.6 


.c 






1.4 


E 


1.2 


*(/) 





600 






CD 




^ 




CD 


400 


D) 




CO 




ak 


200 


o 








m 






e 

4 





M 


= 4 


□ 


M 


= 6 


O 


M 


= 8 



I I 



800 1000 1200 1400 1600 16 
K 



2000 2200 2400 2600 



1 2 3 4 5 6 7 

Breakage location 



where W is a matrix of dimension (M — 1) x (M — 1) 
such that 

Wns = V,,n+i - V;,„ , 1 < n < M, 2 < s < M . (37) 
Here, we have again 

M 



J2 ^ns 



s=2 



A., 



1 



1 < n < Af - 1 



(38) 



For an open-ended string, Eqs (I25p and (pS)) in 

the previous section are modified to: 



(fc) 


= V2Kb 














— rv 



M 
M 

T 



(39) 
(40) 
(41) 



The demonstrations of the above equalities are very sim- 
ilar to those presented in Appendix lA 3[ lA 41 and lA 51 are 
therefore omitted. 

Employing Eq. (|2ip , we arrive at the following asymp- 
totic (as K — > cx)) prediction for the MFBT 



where 



A = 



••• -1 1 



(33) 



with g remains the same as in Eq. ([2]). 

With VDV^ = A, the diagonal elements of D are 



7r(s- 1) 
M 



and V is now defined by 



Vns = C COS 



7r(n- l/2)(s- 1) 
M 



(34) 



(35) 



In the above equation, C — y/2/{M + 1) is again the 
normalization factor as in the case of fixed-ended string. 
The proof for this statement is similar to that presented 
in Appendix I A 1 1 and is thus omitted. 

Note that the first normal mode, pi, corresponds to 
all beads moving in unison and is thus of no interest in 
terms of breakage events. We will therefore omit this 
mode in subsequent discussion. In other words, we will 
only consider the set of normal modes {ps '■ 2 < s < M}. 

As in the previous section, we are interested in the 
extension/contraction given by: 



e = Wp 



(36) 



c 



K3/2fo(M - 1) 



(42) 



This prediction is in good agreement with simulations 
as shown in Fig. 01 Note also that according to this 
calculation, the breakage rate of an open-ended string is 
the same as that of a fixed-ended string when AI ^ 1. 

Since c/)^ and |k£'Q^'^^| are identical for all k, Eq. (^5]) 
predicts that all the springs are broken with equal fre- 
quency. In other words, contrary to the case of an fixed- 
ended string, we would expect a flat distribution of break- 
age frequencies across the string, which is indeed shown 
to be the case by simulations (c.f. Fig. 0]). This re- 
sult may be expected as, unlike in the case of a string 
with fixed ends, the extremal springs here are again con- 
necting two fluctuating beads. We therefore expect that 
their breakage frequencies would be similar to those for 
the springs in the middle of the chain. 

In summary, we have investigated analytically the 
breakage rate for an ID string under thermal fluctuation 
in the heavily damped regime. Our approach is based 
on the theory of the multi-dimensional Kramers escape 
problem, and we have supported our analytical predic- 
tions with numerical simulations. 
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APPENDIX A: PROOFS OF VARIOUS 
IDENTITIES 

1. Claim: = 2/(M + l) 

By definition: 

= ysin^-^ (Al) 

S 

1^/ 2.ns- 



have 



M +1 



ns{n — 1) 
A/ + 1 



(A7) 



2C2^si 



■Ksn , Trs(n — 1) 



sm — sm ■ ^ 

ill + 1 Af + 1 



7rs(2n — 1) TTS 
cos rr : COS ■ 



M+1 



M+l 



M 
s—1 ^ 

M 



(A2) 



M + l 1 27ms , . , 

oE'^^^T^ (A3) 



4. Computations for 



2 2 



Af + 1 



where the second term in Eq. (jA3P is zero as the negative 
terms cancel the positive terms exactly. 



2. Claim: Y.sxT^'^ 
For n = 1, we have 



By Eq. 



M 

A/ 

= (2..)*^ n 

s=l 
M 



1 — cos ■ 



Af + 1 



E 



sin2(7rs/(Af + 1)) 



(4k 



"aT 2 



^ [1 - cos(7rs/(Af + 1)] 



(A4) 



nsin — 
, 2{M+l) 



k'" (Af + 1) 



(A8) 
(A9) 

(AlO) 
(All) 



+ cos(7rs/ (Af + 1)] (A5) where the last equality follows from the identity in \\1 



M 



Af + 1 



(A6) 



The case for n = M follows similarly. 



If n is not 1 or Af , is by definition: 



C2 ^ [sin(7rsn/(Af + 1)) - sin(7rs(n - l)/(Af + 1))]" 



■E 
^'E 

s 

^'E 



1 - cos(7rs/(Af + 1) 
cos2(7rs(2n - l)/(2Af + 2)) sin2(7rs/(2Af + 2)) 



5. Computations for <j)f. 

By rotating the basis, we have a new coordinates, 
y — (j/i, . . . , um), such that the direction (0, . . . , 0, 1) 
corresponds to the direction that is normal to the exit 
boundary at the k-ih. exit point. Let us denote the or- 
thogonal matrix that transforms from this new set of ba- 
sis back to the p-space by R. In particular, we have 



sin2(7r.s/(2Af + 2) 

TTs{2n - 1) 
2Af + 2 



RsM — 



Wks 



(A12) 



Af TTs(2n - 1) 



Af + 1 2 ^ Af + 1 

s 

as the second term in the last equality is zero. 



3. Claim: W^^s equals 1 for n = 1, Af + 1 and 
equals 2 otherwise 

Firstly we prove that 1 f^'' " = 1, Af + 1 and 

is 2 otherwise. The fact that J]^ W^^ = 1 for n = 1, Af-|-1 
follows immediately from Eqs [TlJ For n 7^ 1 , Af -I- 1 , We 



Note that one way to obtain the rest of the matrix ele- 
ments in R is to employ the Gram-Schmidt process. 
In terms of this new coordinates, the potential energy 

is: 



uiy) = ^P^Dp 



= -y^R^DRy 



(A13) 
(A14) 

Restricting to the first Af — 1 dimensions, (j)). is defined 



= det 



([R^DR] 



l<r,s<A/-l 



(A15) 
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We have calculated this quantity for M = 3, 4, . . . , 50 
numerically and we find that the following formula is ex- 
actly satisfied (up to machine rounding errors): 

_ r K*^-iM , /c = l,2,2Af + l,2M + 2 
^fc - I kM-1m/2 , otherwise 

(A16) 

Although we are not able to prove the above equality 
mathematically, we believe that it holds true for all M . 

APPENDIX B: DETAILS OF NUMERICAL 
SIMULATIONS 

Simulations performed for the dynamics of strings with 
fixed ends (free ends) are based on Eq. ([T|) with the ma- 



trix A defined in Eq. ^ (Eq. ([33])). Namely, the posi- 
tions of the beads are updated according to the following 
scheme: 



x(< + At) = x(i) - -Ax(i) At + y —2^ z(t) (Bl) 

where z(t) is a vector with entries given by random num- 
bers drawn from the normal distribution with zero mean 
and a standard deviation of one. The simulations always 
start at the minimal energy configurations and are termi- 
nated when one of the springs' lengths become more than 
1.1 or smaller than 0.9. The parameters in the simula- 
tions are: C, — 10, fc^T = 1, and At is set to be 2 x 10~^. 
For each set of parameters, 1000 runs are performed. 
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